The diffusive instability of kaon condensate in neutron star matter 
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The beta equilibrated dense matter with kaon condensate is analysed with respect to extended 
stability conditions including charge fluctuations. This kind of the diffusive instability, appeared to 
be common property in the kaon condensation case. Results for three different nuclear models are 
presented. 



I. INTRODUCTION 

Following the pioneering work of Kaplan and Nelson 
0, the possibility of having kaon condensation in neu- 
tron stars has attracted a lot of attention during past 
few years. We can divide the various proposed models 
(for a review see 0) into two groups: those coming from 
effective chiral theory and those coming from the meson 
exchange model. The former ones, based on chiral sym- 
metry SU{'i)L X SU{'i)ii, treat kaons as members of the 
Goldstone boson octet and couples them to baryons ac- 
cording to the prescription appropriate for the nonlinear 
representation of the chiral symmetry group . For the 
second group, kaons are not coupled to baryons but only 
to heavy mesons l^i£y'^) ^-nd these in turn inter- 
act with the baryons 0, |^ . Although both approaches 
give similar predictions for the onset of kaon condensa- 
tion, the matter composition remains uncertain in the 
chiral approach due to the shape of the nuclear symme- 
try energy (an important quantity describing the pure 
nucleon interaction in the isovector channel) not being 
well-established. Interesting results, among them the 
prediction of the complete protonization of the stellar 
core, were presented in |6|. However, there is also a qual- 
itative difference between the chiral and meson exchange 
models from a thermodynamical point of view which we 
would like to stress in the present paper. 

Since kaon-nucleon interactions are strongly attractive, 
the presence of a kaon condensate in dense matter makes 
the EOS softer. Because the K-N coupling strength is not 
well known, usually a range of values for it are taken into 
account. For a sufficiently strong, but still reasonable, 
K-N coupling the pressure-energy relation, p{£), exhibits 
a region where the compressibility of the matter is neg- 
ative. This is the signal of phase separation and means 
that one needs to make a kind of " construction" which en- 
sures that the compressibility is positive or at least zero 
(as in the Maxwell construction for the van der Walls 
equation of state). For neutron star matter, the correct 
approach consists of applying the so-called Gibbs con- 



struction. The system may separate into different phases 
provided that their pressures and all relevant chemical 
potentials are equal (the Gibbs conditions). The relevant 
chemical potentials are those connected with quantum 
numbers which are conserved by the interactions main- 
taining the equilibrium. In neutron star matter there 
are two such quantities: the baryon number B and the 
charge Q. This means that the corresponding chemical 
potentials, for neutrons /x„ and for electrons /Xe, must be 
continuous across the phase boundary 



(1) 



These equalities state that the phases are in chemical 
equilibrium - there is no flow of baryon number or charge 
across the phase boundary. Similarly, pressure equality 



pi ^ pii 



(2) 
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states the condition for mechanical equilibrium - the 
phase boundary stays at rest with respect to the matter. 
If the Gibbs conditions \iV2\ can be fulfilled, the system 
forms a mixture of two phases with different densities 
of charge and baryon number, occupying the volume in 
different proportions. This scenario, in the kaon conden- 
sation case, was successfully presented in 0, where the 
meson exchange model was used. Mixed phase formation 
has never been demonstrated for chiral models, although 
the possibility of phase separation was noticed in 
where the Maxwell construction was used to remove the 
negative compressibility region in the EOS. In the case of 
neutron star matter, the Maxwell construction is usually 
made under the assumption that the matter is locally 
neutral. This means that the construction ensures only 
the equality of P and /x„ whereas is different for the 
different phases, so that the EOS obtained is not stable. 
The primary objective of this work was to carry out a 
correct construction to enforce the Gibbs conditions in 
the chiral model and get a stable EOS. However, a thor- 
ough analysis shows that kaon condensation in the chiral 
approach leads to an EOS which is completely unstable, 
independent of the strength of the kaon-nucleon coupling. 
The same kind of stability considerations were also ap- 
plied to the meson exchange approach and we found that 
the stability of an EOS with a kaon condensate depends 
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not on the strength of K-N interactions but on the man- 
ner in which they are reahzed. 

This paper is organized as follows: in Sec^l the more 
general stability conditions for beta stable matter are for- 
mulated. In Sec lIIII it is shown that it is not possible to 
make the Gibbs construction for the chiral model and 
that the stability conditions are not fulfilled. Then, in 
Sec0 the same analysis is repeated for two represen- 
tative meson exchange models, one for which the Gibbs 
construction had been performed earlier by other authors 
and one for which that was not possible. In the final sec- 
tion, we discuss the results obtained and the properties 
of Lagrangians which allow the instability to be avoided. 



II. INTRINSIC STABILITY OF A SINGLE 
PHASE UNDER BETA EQUILIBRIUM 

The basic question in all neutron star matter research 
is the shape of the relationship between pressure and en- 
ergy density p — p{e), usually called the Equation of 
State. At zero temperature, the state of neutron star 
matter should be uniquely described by quantities which 
are conserved by processes leading to the equilibrium. 
For matter above nuclear density uq = 0.16 fm'^ we have 
the nucleon beta cycles 



p + I + i^i p + I 



where I 



e, /i 



which are governed by weak interactions. The baryon 
number B is conserved by this type of process so that 
the energy density e and pressure P should be functions 
only of the baryon number density ub = B/V. How- 
ever, this common statement includes the implicit as- 
sumption that the matter is electrically neutral and spa- 
tially homogeneous. The star as a whole must obviously 
be electrically neutral but the matter does not need to 
be locally neutral. There is no reason why the matter 
needs to be homogeneous - it may separate into two or 
more phases with different charge densities in such a way 
that the total volume on a large scale is neutral. This 
scenario was presented in detail in The second con- 
served quantity, electric charge, should then be treated 
on the same basis as baryon number or, in other words, 
the thermodynamic state of a given phase is described 
not by one but two quantities: baryon number B and 
charge Q, where Q is the sum over all charge carriers 0| 
Q = —Np + Ne + Nf^ + . . . Starting from this point, we will 
formulate correct stability conditions for a system which 
is not neutral, in general, but may have local charge den- 
sity different from 0. Consider a given phase with vol- 
ume V: its energy U is the function of V, B and Q, i.e. 
U — U{V, B, Q). To investigate the intrinsic stability of 
single phase it is more convenient to introduce intensive 
rather than extensive quantities: 

u = U/B, v^V/B, q^Q/B 



and then, the energy per baryon becomes a function of 
two variables 

u = u(v, q) 

and the first law of thermodynamics takes the form 

du = —Pdv + n dq , (3) 
where the pressure and electrical chemical potential are: 

' du 



P : 



du 
dv 



(4) 



From the principle of minimum energy, it can be deduced 
that the phase is intrinsically stable (i.e. it does not 
separate into different phases) if and only if the energy 
per baryon is a convex function of its variables v and q (for 
a detailed analysis, see 0). In terms of its differentials 

d^u = ^Uyy dv^ + Uyq dv dq + \uqq dq^ > 0. (5) 

The sufficient condition for the quadratic form to be pos- 
itive is that its determinant must be positive and at least 
one of derivatives Uy^ or Uqq must be positive: 



Uyy Uyq 



> and ( Uyy > or Uqq > 0) (6) 



Thermodynamical identities applied to the above second 
derivatives allow us to express the positivity condition in 
terms of two equivalent pairs of inequalities: 



dP 
dv 



/dP 



>0 ^ 



du 
dq 



or 



>0 



d^i 
dq 



> 



> 



(7) 



(8) 



The pressure derivatives correspond to the well-known 
compressibilities of matter, one with constant charge and 
the other with constant chemical potential: 



^ "Kl^ ) «M = -7 ( ^ ) (9) 



1 / dv 
V \'dP 



The positivity of the compressibility, usually referred to 
as mechanical stability, here means that baryon den- 
sity fluctuations are stable. The physical content of 
the fj, variation in Q and ^ becomes clear if we look 
carefully at their dimensions. As fi has dimensions 
of energy/charge, these derivatives have dimensions of 
energy /charge^ which is the inverse of the unit of elec- 
trical capacitance. It is useful to introduce the electrical 
capacitance of matter: x- Obviously, we must specify un- 
der which constraints it is taken - under constant pressure 
or constant volume: 



dp 



dq 



(10) 
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When X is positive, it means the system is diffusively sta- 
ble, i.e. the charge fluctuations are stable. If x changes 
its sign, the charge fluctuations become infinite which 
means that the phase can be no longer homogeneous and 
separates itself into parts with different charge density. 
So, finally, it may be concluded that for stable matter, 
the compressibility and electrical capacitance must be 
positive, what may be expressed in two ways: 

f '^g > > 

I Xp > I X, > ^^^^ 

We can notice here some similarities to the one com- 
ponent system with non-zero temperature, commonly 
treated in textbooks on thermodynamics. For this kind 
of system, the internal energy [/ is a function of volume, 
entropy and particle number: U {V, S, N). If we make the 
following replacements: 

q ^ s 
^ T 

i.e. the charge per baryon plays the role of specific en- 
tropy s = S/N and the chemical potential plays the role 
of temperature, then the electrical capacitance Xi corre- 
sponds to the specific heat: 

Instead of compressibility under fixed q or fi, we here 
have compressibility under constant T or s and, in the 
same way, convexity of the internal energy, which is now 
a function of specific volume and entropy u{v, s), requires 
the standard compressibility and specific heat to be pos- 
itive: 

J Kj, > ( Ks > 

{ Cy > \Cp > 

which are well known sufficient conditions for stability of 
a one component system at non-zero temperature. 

In the next section we stress that the positivity require- 
ments on the electric capacitance of matter may not al- 
ways be ensured. This property is strictly connected with 
the interactions used for the description of the dense nu- 
clear matter. 



III. THE CHIRAL MODEL 

To investigate phase stability in the chiral approach, 
we focus on the model used in where kaon-nucleon 
interactions are described by the Lagrangian coming from 
chiral effective field theory 

= i^Trd^.Ud'^U'' + TrB{ij''D^^mB)B (12) 
+^^TrB7^75[A., B] + OTrB7^75{^^, B} 
+cTrM{U + C/t) + aiTrB{(,M(. + ^''M^^)B 
+a2TrBB{^M^ + ^''M£.^) 
+as.TrBB Tr{^M^ + S}M£}). 



The only parameter in which is not well determined 
is as, which is connected with the kaon-nucleon sigma 
term Y^kn — \ms{N\uu + ss\N) — —{\ai + a2 + 
2a3)TOs- Reasonable values for a^ma are in the range 
—310 MeV < a^ms < —134 MeV which corresponds to 
a strangeness content for the proton within the range 
0.2 > {ss)p/{uu + dd + ss)p > lll|. The value of this 
parameter does not appear to be of great importance for 
the following analysis of matter stability. The energy 
density coming from is found by means of Baym's 
theorem, which states that {K~) — ■^exp(— /^^t) (as- 
suming only the s-wave kaon-nucleon interaction): 

£KN = sin^ e + 2{mlf - n^KNix)) sin^ ^ (13) 

where 

Ykn{x) = -{aix + 02 + 2a^)ms. 

The kaon-nucleon contribution must be properly mini- 
mized with respect to 6*, and this leads to the following 
equation for 9 

cos 6* = ^ ^ {'m^Kf - nYjKN{x) - \iiKn{\^-x)) . 

(14) 

For the interactions in the pure nucleonic sector we adopt 
the phenomenological formula ^2] for the energy density 

£NN = ^E'^p^noi^f/'' + V{n) + mn + n{l - 2xfEs{n), 

(15) 

where 

Es{n) = ^(22/3-l)4")((^)2/3-F(;^)) 

+£;.(no)i^(^). (16) 

There are two potential contributions in H15|) com- 
ing from the isoscalar and isovector sectors of the nu- 
cleon interaction, V{rL) and Es{n) respectively. V{n) is 
parametrized to get correct values for the compression 
modulus at the saturation point ~ 240 MeV and it 
mainly determines the stiffness of the EOS. The nuclear 
symmetry energy Eg (n) at the saturation point takes the 
value 30 MeV. The potential contribution F{u) is an 
unknown function which determines the behaviour of E^ 
at higher densities. The importance of the symmetry 
energy for the chemical composition of dense matter in 
the kaon condensate case was presented in Here we 
put F{u) = u to mimic typical results from relativistic 
mean field theory (RMF) since we are going to compare 
the chiral model to the RMF approach only as regards 
the difference in the kaon-nucleon interaction. The total 
energy density is a sum 

£ = £KN + £NN + ^lep, (17) 

where the leptonic contribution is expressed in the stan- 
dard way by means of chemical potentials for electrons 
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and muons 



£iep = £^ + ^t9i\/vitJ-l -ml)/m^) , (18) 
where the function g[t) is: 

g{t) = -^((2^3 + t)\/l + t^ ~ arsinhi) , (19) 



77(x) 



a; ; a; > 
; a; < 



From the total energy of the system we are able to find 
all of the necessary thermodynamic quantities such as the 
nucleon chemical potentials and the pressure 

= Ef{^y/^ + m + V'{n)+n{l-2xfE',{n) (20) 



+ {l-Ax'^)Es{n) - (2I]A'„ + MK)sin2 2' 
= E^P{^y'' + m + V'{n)+n{l-2xfE',{n) (21) 

Q 

-{l-2x){i-2x)Es{n) - (2Eifp + 2/Ltif)sin2-, 
P = -e + ^^fiiUi = 



+ nV'{n) + n^{l - 2xyE'^{n) 
+ l-f^ilsm^e-2mlfsm^^ 



(22) 



+ 



where for the chemical potentials it is useful to intro- 
duce kaon-neutron and kaon-proton sigma terms: T.Kn = 

^Kn{0),'^Kp — S/f7v(l)- 

Now, the total energy and pressure are functions of five 
variables (the number densities nn,np, fj.K, iJ^e, ^J■fl cor- 
responding to the five species of particles p,n, K~ ,e, ^ 
present in the system). By the use of the beta equilib- 
rium conditions 



(23) 



we may reduce these to two independent variables as 
mentioned in section ^] The above equations also allow 
us to show that 



Me 



du 



which justifies identifying /ie with the charge chemical 
potential appearing in |0J 

M = Me- 

Which two of the five variables are chosen is a matter of 
convenience. For instance, in order to make the Gibbs 
construction we need the pressure as function of the 



chemical potentials P(/i„,/i) for finding a state preserv- 
ing the Gibbs conditions (|1I2I) for the two-phase coexis- 
tence. Usually, the compressibility of the neutral matter 
is first tested (§^)^_q > 0, and following this procedure, 
for the full stability considerations we have chosen the 
pair of inequalities expressed by O . For the ^ derivative 
in |7J it is useful to express the pressure as a function of 
/X and q, where 



nx + f fisin 6 + n{l + x) sin % + ni 



ep 



(24) 



and then to look at the pressure contours in the q-^ plot. 
Using the condition of beta equilibrium (|23|l , the pressure 
in Ea. H22|) can be written as a function of n, x, 9, fi. These 
quantities were found numerically for given P and q by 
solving Eqs H14I22I24|I . with the proton fraction under 
beta equilibrium being 



1 /i(l + cos 0) + 2aims(l — cos 6*) 

2 lGEs{n) ' 



(25) 



IV. PHASE STABILITY IN CHIRAL MODEL 

For the kaon condensation presented in the 
pressure-density relation revealed a region of negative 
compressibility when the matter was locally neutral 
q = 0. This unstable region was removed by means of 
the Maxwell construction which, as mentioned in the in- 
troduction, does not preserve continuity of the charge 
chemical potential. By abandoning the local neutrality 
of matter one may to try to perform the Gibbs construc- 
tion as presented in Fig^ where the pressure is shown as 
a function of the baryon density and the charge chemical 
potential and also /i discontinuity is shown in the case of 
the Maxwell construction. From this figure one can see 
that the pressure of the kaon phase is always lower than 
that of the normal phase. The only points where the two 
pressures are equal is on the 9 = line, where the con- 
densate starts to grow. However, at this line these two 
surfaces are tangent, or the pressure gradients for both 
phases are equal: 



N 



dP 



dP 



K 



(26) 



These pressure gradients are just the charge densities, 
according to the thermodynamical relation 



dP = n dfin + nq d/i 



(27) 



where riq = Q/V, so there is no difference in charge den- 
sity between the normal and kaon phases - the transition 
from the normal phase to the kaon phase is of second 
order. The lack of charge difference between the phases 
makes it impossible to have formation of the mixed phase 
along 9 = 0. There is no other region in the (/i„, /i) plane 
where the pressures for the two phases might be equal. 
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FIG. 1: Left panel: the pressure sheets for the normal phase (upper) and the kaon phase (lower) for a^ma = —230 MeV. The 
thin dashed line corresponds to 6 = where condensation starts and the two sheets are tangent. The thick line represents the 
neutral matter where its dashed part is the projection of the line which, in fact, lies on the lower sheet corresponding to the 
kaon phase. The right-hand panel shows the Maxwell construction which preserves the continuity only of P and 



This is a signal of inconsistency between pressure and 
energy for the kaon phase. Thermodynamic principles 
state that the system prefers the phase having greater 
pressure for given chemical potentials ("pressure maxi- 
mization") and simultaneously smaller energy for given 
particle numbers (" energy minimization" ) . Going along 
the neutrality line with increasing density, after the un- 
stable region is passed, the line should have entered into 
region where the mixed phase disappears and the pure 
kaon phase is again preferable. Although for all den- 
sities the kaon phase has lower energy, suggesting that 
it is the preferred phase, its pressure never exceeds 
. The fundamental reason for this is that the kaon 
phase does not represent a stable system at all. This 
conclusion is confirmed by testing the second derivative 
of the energy, Xp^ which is shown by means of constant 
pressure contours in the q-^ plot, in Fig|21 The isobars 
are represented keeping in mind the physical constraints: 
< X < 1 and — 1 < cosO < 1, and so usually they do not 
cover the whole area of the plot. As one can see, the iso- 
bar slope for the kaon phase is negative for the neutral 
phase q = 0. Moreover, a negative Xp is also observed 
for almost all values of q. For comparison, the isobars for 
the normal phase are also plotted and they always ex- 
hibit a positive value of Xp which means that the normal 
phase is always stable. Such behaviour is not an effect of 
the kaon-nucleon coupling being too strong. A negative 
Xp is also observed for the minimal kaon-nucleon cou- 
pling corresponding to the vanishing strangeness content 
of proton where 037715 — —134 MeV. The plots for this 
minimal coupling are presented in Fig|31 The pressure- 
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FIG. 2: Isobars on the g-^ plot for the normal phase (dashed 
lines) and the kaon phase (solid lines) . The values of pressure 
indicated are expressed in MeV/fm"^. 



density relation for the neutral matter does not reveal 
any region of negative compressibility, as seen in the left 
panel in Fig|3| Nevertheless the electrical capacity Xp 
is still negative for g = 0. This means that the matter, 
although being mechanically stable, does not fulfil the 
diffusive stability condition (TJ. Also, for most values of 
(7, kaonic matter is diffusively unstable except in a small 
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region for low values of pressure and only positive values 
of q. At this point we conclude that, within the frame- 
work of the chiral model, matter with a kaon condensate 
is intrinsically unstable with respect to charge fluctua- 
tions for any values of K-N coupling. 



Apart from the usual Weinberg- Tomozawa term (the sec- 
ond term in £^^^), the covariant derivative introduces 
additional terms which are quadratic in the meson fields 
and Rf^. These additional terms modify the equation 
of motion for the vector mesons in such a way that it 
couples to the conserved current 



V. STABILITY IN THE MESON EXCHANGE 
MODEL 

Having discussed the stability considerations for the 
chiral model we now want to address the same question 
for the second class of models used for treating kaon con- 
densation - the meson exchange models. The common 
feature of this class of models is that kaons are coupled 
to vector or scalar mesons but not directly to baryons. 
The way of making this coupling is not unique. Let us fo- 
cus on two different ways of coupling presented by Knor- 
ren, Prakash and Ellis in Q and by Glendenning and 
Schaffner-Bielich in (2| . In the following, we refer to them 
as the KPE and GS approaches. For the GS model, the 
Gibbs construction was successfully carried out whereas 
for KPE it was not possible to do this as was later men- 
tioned in The two approaches differ regarding the 
pure nucleon interactions as well as the kaon-nucleon in- 
teractions. In order to underline the role played by the 
difference in the kaon-nucleon sector we have formulated 
two models which are identical in their nucleon part but 
different in the kaon part. Hence, for both of them we 
keep the same Lagrangian for the non-strange sector 

+ \d^ad^a - \m'y - U{a) 



(28) 



where N denotes nucleons and cr, V, R are the scalar, vec- 
tor and isovector meson fields (for details see i^] ) . In the 
kaonic sector we used two different Lagrangians: 



rGS 



(29) 



(30) 



where the covariant derivative is 

'Df,=dfj, + ig^jKVf, + igpKTKR^, (31) 
The two Lagrangians have different mass terms: 



m 



<GS) 
K 



m 



'■{KPE) 
K 



mK - gaKcr 



^x ~ gcrKfriKCF 



(32) 



(33) 



d^V^" + m^V^ = g^NN-i^'N - g^KJ'ti (34) 

where the kaon current Jj^, coming from the symmetry 
exp±*"i^±, is 

= iK-d^K+ -idpK-K+ 

-2g^KV^K-K+ - 2gpKTKRt.K-K+ (35) 

In this scheme, the divergence of the vector field van- 
ishes and so it is possible to fulfil the Ward identity for 
a vector field in a medium, as pointed out by Schaffner 
and Mishustin in In the KPE Lagrangian, the cou- 
pling to the scalar a is only linear whereas for the GS 
Lagrangian a quadratic term in a is also present. 

For both models presented here, we use the same 
parametrization for the nucleon part (|28|l which gives the 
properties of the saturation point as follows: the bind- 
ing energy at uq = 0.153 fm^^ is E/A — —16.3 MeV, 
the symmetry energy Es{nQ) = 32.5 MeV, the compress- 
ibility Kq = 240 MeV, and the nucleon effective mass 
m*/m = 0.78. As shown in 0, the density dependence 
of the symmetry energy in the RMF models is domi- 
nated by a linear increase with density which is in agree- 
ment with the Es{n) introduced in the chiral approach 
(|16ll . In order to keep the same strength for the K-N in- 
teractions in both models, one may compare the optical 
potential in symmetric nuclear matter (i?o,3 = 0) coming 
from (1201) and 123 



ttGS 



-gaKm„ - g^KVo + 



iguKVo)'^ + igaK(T) 



2mK 



U§^^ = --gcrK-rria - gujKVo 



(36) 
(37) 



At nuclear density no, where Uk is fixed, the quadratic 
terms in U'~^^ contribute no more than a few percent to 
the total and so we neglect them, as was also done in 9] . 
The only difference is then in the scalar coupling constant 
ga-K which is twice as large in the KPE parametrization 
as in the GS one. The value of the optical potential may 
also be used to relate the K-N coupling strength to the 
chiral model for which Uk is 



jjX _ ^KNK^j 



3n 



(38) 



A strangeness content of the proton within the range 0- 
0.2 then corresponds to Uk being between -73 and -120 
MeV. To keep contact with the results of Glendenning 
and Schaffner, we use the values -80 and -120 MeV in 
our further analysis. 
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FIG. 3: Kaon phase instability for minimal K-N coupling. 



The models formulated above are treated in the spirit 
of RMF theory, taking into account that the kaon mean 
field value is — ^ exp^'^* . In order to find the pres- 
sure contours in the q-^ plot, one has to solve the relevant 
equation of motion including the beta stable matter (see 
[3) for given P and charge per baryon q, where 

q = 39 

n 

and the kaon number density is 

HK = 2{fi + g^KVo + gpKRo,z)K-K+ (40) 

First, we present the stability considerations for the GS 
model. In Fig0|one may compare the results for the two 
values of the optical potential. For the Uk — —80 MeV 
case, when the neutral kaon phase appears its compress- 
ibility is positive at all densities (the left-hand plot). The 
kaon phase is also diffusively stable, its isobars having a 
positive slope so that that Xp > 0- For the Uk = —120 
MeV case, just after the onset of condensation the com- 
pressibility of the neutral phase and Xp ^^re both nega- 
tive. To be more explicit, the isobar P = 89 MeV/fm'^, 
corresponding the critical density Uc is plotted to show 
that it keeps the negative slope for q = 0. The diffu- 
sive instability means that the charge fluctuations in- 
crease infinitely leading to phase separation. The charge 
is unequally deposited into the two different phases which 
can coexist only if the conditions of mutual stability, ex- 
pressed by equations and arc fulfilled. The coex- 
istence line (absent in the chiral approach) can be found 
for the GS model from the equation 

P^(Ai„,/^) = P^(Ai„,/i) • (41) 

In the q-^ plot, the coexistence line corresponds to two 
lines: N1-N2 and K1-K2 because they are located on two 



different pressure sheets P^ {q, ^) and P^{q,^). Pro- 
ceeding with increasing baryon density, at the point Nl 
the system splits up into a nucleon phase and a kaon 
phase (the point Kl). At the point Kl, the kaon phase 
is maximally charged with a negative charge (q is posi- 
tive for negative charge). Increasing the density further 
from Kl to K2, the kaon phase loses its charge in favour 
of the nucleon phase, becoming neutral at the point K2. 
At this point the pure kaon phase starts to be preferred, 
the nucleon phase vanishes and again the system is in 
a one-phase state. We want to point out that the lines 
for both the normal and (particularly) the kaon phase 
pass through a region where the phases are intrinsically 
stable, with the corresponding isobar slopes always be- 
ing positive, ensuring that Xp > 0. In the mixed phase 
region, the charge q and the baryon number density n 
are unequally distributed on the micro-scale. Within a 
sufficiently large volume V , one may observe droplets of 
one phase immersed in the other. The volume fraction 
a = /V for the two phases is found from the require- 
ment of the global neutrality: 

+ = « g^n^ + (1 - a) g V = 0. (42) 

We may then define the average baryon number density 
for the mixed phase region 

n^an^ + {l~a)n^ (43) 

and represent the other thermodynamical quantities 
{P,fj,,e etc.) as functions of the one variable n (see the 
dotted lines in the left-hand panel in Fig^. The above 
analysis confirms that for the GS model the intrinsic sta- 
bility conditions for a given phase Q and the mutual 
stability conditions for the mixed phase H1I2|) are simul- 
taneously fulfilled. As we will see in the following, this 
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FIG. 4: Kaon phase stability for the GS modeL The left-hand panels show quantities for neutral phase, the dotted lines 
representing the values for the mixed phase region. In the right-hand panels, isobars for the normal and kaon phases are 
plotted. For the Uk = —120 case, the phase separation is also shown. 



does not happen for other meson exchange models like 
the KPE one. The results for the KPE model are very 
similar to those for the chiral model. Results for the two 
cases Uk— — —120, —90 are shown in Fig[S| In the first 
case, negative compressibility appears for neutral phase 
just after the critical density. It is impossible to find a 
region in the q-fi plane where the kaon phase would be 
diffusively stable. Even at higher densities, correspond- 
ing to negative values of /x, when the compressibility of 
q = phase again becomes positive the kaon phase is still 
diffusively unstable. As in the chiral case, this property 
of the KPE model remains unchanged also for smaller 
values of the kaon-nucleon coupling. As one may see in 
FigEl for Uk — —90 (we took this value instead of -80 
so that the effect is more visible), in spite of positive Kq 
for the whole range of densities considered, the electrical 



capacitance Xp is always negative apart from in a small 
range of q around 0.6 which corresponds very small values 
of the pressure. 

VI. SUMMARY AND DISCUSSION 

In this paper we have focused on the diffusive stability 
of kaon condensates within the framework of three dif- 
ferent models of strong interactions: one chiral and two 
meson-exchange based models (here referred to as GS 
and KPE). From this analysis, one can understand why 
the Gibbs construction, which can be made in the GS 
model, is not applicable to the other models presented 
here. From a thermodynamical point of view, the chiral 
and KPE Lagrangians lead to matter which does not ful- 
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FIG. 5: The kaon phase instabihty for the KPE model with optical potential Uk ~ —120 and —90 MeV. 



fil one of the fundamental conditions for intrinsic phase 
stabihty. This happens independently of the value of the 
kaon-nucleon coupling strength and must be an intrinsic 
property of the model. At this point it is difficult to say 
definitely what are the requirements which should be im- 
posed on a Lagrangian in order to get complete stability. 
In the work by Pons et. al the authors considered 
several meson exchange models and found that making 
the Gibbs construction is possible for some of them. One 
may suppose that the existence of a mixed phase exis- 
tence means that the kaon phase is intrinsically stable, 
but in principle, this should be checked independently for 
the models presented. Looking at the properties of the 
Lagrangian for the admissible models, one may suppose 
that nature of the kaon to vector-meson coupling is the 
most relevant property. Making the Gibbs construction 
was possible only for models with the covariant derivative 
for the kaons given by ea. H31|l . which leads to vanishing 



divergence of the vector field. However this condition 
is not sufficient; as was shown in p^ . there was also a 
model with this property for which the Gibbs construc- 
tion was not possible. It seems that the kaon to scalar 
meson coupling is also important. To reach a more solid 
conclusion on this will require further analysis. 
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